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Abstract We theoretically investigate magnetic properties of a unitary Fermi gas in a har¬ 
monic trap. Including strong pairing fluctuations within the framework of an extended T- 
matrix approximation (ETMA), as well as effects of a trap potential within the local density 
approximation (LDA), we calculate the local spin susceptibility x(T,r) above the superfluid 
phase transition temperature T c . We show that the formation of preformed singlet Cooper 
pairs anomalously suppresses X{T,r) in the trap center near T c . We also point out that, in 
the unitarity limit, the spin-gap temperature in a uniform Fermi gas can be evaluated from 
the observation of the spatial variation of %(T, r). Since a real ultracold Fermi gas is always 
in a trap potential, our results would be useful for the study of how this spatial inhomogene¬ 
ity affects thermodynamic properties of an ultracold Fermi gas in the BCS-BEC crossover 
region. 

PACS numbers: 03.75.Hh, 05.30.Fk, 67.85.Lm. 


1 Introduction 

An ultracold Fermi gas provides us the unique opportunity that we can systematically study 
physical properties of a many fermion system at various interaction strengths, by adjusting 
the threshold energy of a Feshbach resonanc e 1 - 2 i 3 i 4 . Indeed, by using this advantage, the so- 
called BCS (Bardeen-Cooper-Shrieffer)-BEC (Bose-Einstein condensation) crossover— 1 ' 7 — 
has experimentally been realized in 40 K— and Fermi gases, where a BCS-type 

Fermi superfluid continuously changes into the BEC of tightly bound molecules, with in¬ 
creasing the strength of a pairing interaction. In this sense, we can now deal with a Fermi 
superfluid and a Bose superfluid in a unified manner. 

Recently, the spin susceptibility has become accessible in the BCS-BEC crossover regime 
of an ultracold Fermi ga s 14 i 15 i 16 . Here, “spin” is actually pseudospin describing one of the 
two atomic hyperfine states contributing to the pair formation. Using this thermodynamic 
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quantity, we can examine to what extent the spin degrees of freedom are active in the BCS- 
BEC crossover region. Theoretically, the possibility of the so-called spin-gap phenomenon 
has been discussed in the crossover region near the superfluid phase transition tempera¬ 
ture T c 17 ] 18 i 19 i 20 i 2 i , where the spin susceptibility is anomalously suppressed by preformed 
spin-singlet Cooper pairs. Since preformed Cooper pairs also cause the pseudogap phe¬ 
nomenon 22 ’—' 24 is4s& (where the single-particle density of states exhibits a gap-like structure 
even in the normal state), the spin-gap phenomenon and pseudogap phenomenon are deeply 
related to each other in the cold Fermi gas system. 

So far, the spin susceptibility has theoretically been discussed in a uniform Fermi gas 18 - 19 ' 20 - 21 
although a real ultracold Fermi gas is always prepared in a trap potential. In this paper, thus, 
taking this realistic situation into account, we study how spatially inhomogeneous pairing 
fluctuations affect the spin-gap phenomenon in a trapped unitary Fermi gas. For this pur¬ 
pose, we employ the extended T -matrix approximation (ETMA) developed in the uniform 
syste m 18 ' 21 , to include effects of a harmonic trap within the local density approximation 
(LDA ) 23 1 26 . In a uniform Fermi gas, it has been shown that ETMA correctly describes the 
BCS-BEC crossover behavior of the spin susceptibility^^!, which makes us expect that 
this strong-coupling theory is also valid for the trapped case. We briefly note that the or¬ 
dinary T-matrix approximation 22 - 25 - 27 , as well as the strong-coupling theory developed by 
Nozieres and Schmitt-Rinfc^2, are known to unphysically give negative spin susceptibility 
in the BCS-BEC crossover region, although these theories have successfully explained vari¬ 
ous many-body phenomena in the BCS-BEC crossover region. Using the combined ETMA 
with LDA, we calculate the local spin susceptibility x{T,r) in the normal state near T c . 
Throughout this paper, we take h = k b = 1, for simplicity. 


2 Formulation 


We consider a two-component Fermi gas, described by the BCS Hamiltonian, 
H = H^ c l,c c v,o-u £ C p+ f,t c - P + f,! c - P '+§,-i c P '+f 

p-tf p-p 


( 1 ) 


where c p a is a creation operator of a Fermi atom with pseudospin a =f, J,. | p a =p 2 / (2m) — 
p — ah is the kinetic energy in the cr-spin component, which is measured from the Fermi 
chemical potential p, where m is an atomic mass, and h is an infinitesimally small effective 
magnetic field to calculate the spin susceptibility. The pairing interaction —U is assumed to 
be tunable. The unitarity limit (which we are dealing with in this paper) is characterized by 
the vanishing inverse s-wave scattering length = 0), which is related to the interaction 
strength —U as 


4 7ia s 


U 


l-t/I'p 1 


Pc m ’ 


( 2 ) 


where p c is a cut-off momentum. 

In LDA, effects of a harmonic trap potential V(r) = mQrr 2 /2 can be conveniently in¬ 
corporated into the theory by simply replacing the chemical potential p a = p + ah with 
the position-dependent one p a (r) = p 0 — V (r) 23 i 26 , where Q is a trap frequency. The LDA 
single-particle thermal Green’s function then has the form. 


G a ip,i(O n ,r) = 


Opty (r) Try f p. io ) n . r ) 


( 3 ) 









Title Suppressed Due to Excessive Length 


3 


(a) 

(b) 





Fig. 1 (Color online) (a) Self-energy correction Z a (p ,i(Q n ,r) in the extended T -matrix approximation 
(ETMA). (b) Particle-particle scattering matrix F(q,z’v n ,r). The solid line and double-solid line represent 
the bare Green’s function and the ETMA Green’s function G a , respectively. The wavy line describes the 
pairing interaction —U. 


where co n is the fermion Matsubara frequency, and £p iC7 (r) — P 2 /{2m) — flair). The LDA 
self-energy ,/co /7 ,r) describes fluctuation corrections to single-particle Fermi excita¬ 
tions. In ETMA, it is diagrammatically described as Fig. |T] which gives, 

E a (p,ico,„r) = T £ r(q,iv n ,r)G- a (q-p,iv n -i(O n ,r). (4) 

q ,iv„ 


Here, v„ is the boson Matsubara frequency, and — <7 means the opposite component to <7- 
spin. We briefly note that the ordinary T-matrix approximation 22 - 25 - 27 is immediately re¬ 
produced by simply replacing the ETMA Green’s function G- a in Eq. 0 with the bare 
one, 


G 0 _ a (pJm n ,r) = — 

1 ( 0 , 


1 

%p-o( r ) 


( 5 ) 


In Eq. 0, r (q, iv„, r) is the particle-particle scattering matrix, given by 


r(q ,iv n ,r) 


-U 

1 -UII(q,iv n ,r)’ 


( 6 ) 


where 

n(q,iv n ,r) = T £ G°(p + q,iC 0 „ +;'v„,r)G°(-p,-/ft)„,r) (7) 

is the lowest-order pair-correlation function, describing fluctuations in the Cooper channel. 

In the present formalism, the superfluid phase transition temperature T c is determined 
from the condition that the Thouless criterion is satisfied at the trap center as F (q = 0, iv„ = 
0,r = 0) _1 = 0 23 i 26 . As usual, we solve this equation, together with the equation for the 
total number N of Fermi atoms, given by 

V = E fclnifj(r) , (8) 

G J 


to self-consistently determine T c and /i. Here, 

n a (r) = T £ G a (p,iCO n ,r) 

P,iO)n 


(9) 
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Fig. 2 (Color online) (a) Calculated local spin susceptibility %(T,r) of a trapped unitary Fermi gas, as a 
function of the spatial position r measured from the trap center. At each temperature, the spatial position 
^sg(T) at which %(r,r) takes a maximal value is shown as the filled circle. The horizontal dotted line shows 
the maximal value (~ 0.42) of the sca led spin susceptibility in the case of a uniform unitary Fermi gas at T c 
shown in panel (b)^i. R f = yj2ep/m/Q is the Thomas-Fermi radius, where £p is the LDA Fermi energy at 
the trap center. X^iT — 0,r) is given in Eq. ced (b) Spin susceptibility Xum{T) in a uniform unitary Fermi 
gas 18 i 21 . The filled circle shows the spin-gap temperature T s ~ 0.377j? at which Xuni{T) takes the maximal 
value, Xuni (T) — 0.42^°, where x° is the spin susceptibility of a free Fermi gas at T = 0. is the Fermi 
temperature in a uniform Fermi gas. 


is the local number density of Fermi atoms with a spin. Above r c , we only solve the LDA 
number equation <[8]), to determine the chemical potential /i. 

The local spin susceptibility x(T, r) is calculated from, 


M ’ ’ dli h->0 h 


( 10 ) 


In this paper, we numerically evaluate Eq. dUB- by taking a small but finite value of h. 


3 Inhomogeneous spin-gap phenomenon in a trapped unitary Fermi gas 

Figure [2f a) shows the local spin susceptibility x(T,r) in a trapped unitary Fermi gas above 
T c . Here, x(T,r) is normalized by the zero-temperature spin susceptibility X°(T = 0,r) in 
an assumed uniform free Fermi gas with the particle density being equal to the density at r 
in the trapped case, given by 

z °(r = 0,r) = 3m(37T 2 ) - 3n(r)3, (11) 

where n(r) = n-f(r) +nj,(r). Since the density profile monotonically decreases as one goes 
away from the trap center (See Fig.[3ja).), pairing fluctuations become weak around the edge 
of the gas cloud even at T c . On the other hand, atoms feel a high scaled-temperature T/Tp(r) 
around the edge of the gas cloud, because the LDA local Fermi temperature, 

T F (r) = [3n 2 n(r)]^ 3 /2m, (12) 

is low in the low-density region. (See Fig. [3jb).) As a result, the local spin susceptibility 
X{T,r) is suppressed thermally around the edge of the gas cloud, as in the case of a simple 
free Fermi gas at high temperatures. Thus, one has x(T.r)/x°{T = 0, r) <§; 1 in this spatial 
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Fig. 3 (Color online) (a) Density profile n(r ) in a trapped ultracold Fermi gas at various temperatures, (b) 
Scaled temperature T/Tp(r), as a function of r. The inset shows T/Tp(r) magnified around the trap center, 
where r^aiT) is the peak position of %(T, r )//t°(T = 0, r) in Fig.[2la). The horizontal dotted line in the inset 
shows the spin-gap temperature T s /Tp — 0.37 in a uniform Fermi gas at the unitarity. 


region, as seen in Fig.[2ja). This ordinary thermal effect becomes weak, as one approaches 
the trap center, because of the decrease of the scaled temperature T /7p(r), as shown in Fig. 
[3}b). As a result, x{T : r)/%°(T = 0, r) increases, as one approaches the trap center from the 
outer region of the gas cloud. 

However, Fig. [2ja) shows that the scaled spin susceptibility %(T = T c ,r)/x°(T = 0,r) 
is suppressed in the vicinity of the trap center, r < ().46Rp (where R\ is the Thomas Fermi 
radius), in spite of the fact that the scaled temperature T/Tp(r) still decreases with decreas¬ 
ing r in this spatial region (because of the monotonic spatial variation of the density profile 
shown in Fig.[3ja)). Thus, this suppression is not due to the simple thermal effect, but is con¬ 
sidered as the spin-gap phenomenon originating from strong pairing fluctuations enhanced 
in the trap center near T c . Indeed, in the spatial region r < j'sg(T’), where /"sg(D is the po¬ 
sition at which x(T-,r)/x°(T = 0 ,r) takes a maximal value, x(T,r)/x°(T = 0, r) is found 
to increase with increasing the temperature. While this temperature dependence is oppo¬ 
site to the case of a uniform free Fermi gas (where the spin susceptibility monotonically 
decreases with an increase of the temperature), it is consistent with the temperature depen¬ 
dence of the spin susceptibility in the spin-gap regime (T < T s ) of a uniform Fermi ga s— 2 L. 
(See Fig. |2jb).) As shown in Fig.[2ja), the spatial region, r < rsdT), becomes narrower at 
higher temperatures, to eventually vanish at T ~ 1.33r c , reflecting the weakening of pairing 
fluctuations. 

A uniform Fermi gas at the unitarity is known to exhibit the so-called universal thermo¬ 
dynamics 31 - 32 - 3 -\ where the Fermi energy £p (or equivalently the Fermi temperature 7p) is 
the unique energy scale, because of the vanishing inverse scattering length a ~ 1 =0. In the 
present trapped case, the scaled local spin susceptibility in LDA is expected to behave as, 

-0^L- )=X {T/T ¥ {r)). (13) 

The same universal function X(x) in Eq. ( 113b is also expected in the uniform case as 

Xuni(r) 

x° 


X(T/T P ), 


(14) 
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where Xuni is the spin susceptibility in a uniform unitary Fermi gas, and x° is the zero- 
temperature susceptibility in a uniform free Fermi gas. 7p is the Fermi temperature in a 
uniform free Fermi gas. Using the relation between Eqs. d 1 3t and < 114l ). together with the 
fact that the scaled temperature T/Tp(r) is related to the spatial position through Eq. ( 1 1 2b . 
we can relate the spatial variation of x{T,r)/x°{T = 0,r) in Fig. [3a) to the temperature 
dependence of Xum{T)/x° i n Fig. [2jb). Indeed, the maximal value Xam/X° — 0-42 at the 
spin gap temperature T s /Tp ~ 0.37 in a uniform unitary Fermi gas (Fig. (3b)) just equals 
the peak value of x(T,r)/x°( r > T = 0) at r = rsc(T) in the trapped case (Fig. [3b)), and the 
latter result is independent of the value of T. In addition, the inset in Fig. [3jb) shows that 
the local scaled temperature T /7p(r = rsc(T)) in the trapped case always equals the spin 
gap temperature T s /Tp — 0.37 obtained in the uniform case. These universal results indicate 
that the observations of the spatial variation of the spin susceptibility x(T,r), as well as the 
density profile n(r), in a trapped Fermi gas at the unitarity enable us to evaluate the spin-gap 
temperature T s in a uniform unitary Fermi gas. 

In this regard, we briefly note that the relation between a uniform Fermi gas and a 
trapped one become complicated when a~ s 1 / 0. In this case, the LDA spin susceptibil¬ 
ity x(T- r ) in a trap also depends on (pp(r)a s ) _1 in addition to T/Tp(r), where pp(r) = 
[3W 2 n(r)] 1 / 3 is the LDA local Fermi momentum. As a result, x(T, r) is related to the spin 
susceptibility in a uniform Fermi gas, not only at various scaled temperatures T /7p, but also 
at various interaction strengths (ppa s )~ l , where pp is the Fermi momentum in a uniform 
Fermi gas. 


4 Summary 

To summarize, we have discussed magnetic properties of a unitary Fermi gas in a harmonic 
potential above T c . Including strong pairing fluctuations within the framework of the ex¬ 
tended T -matrix approximation (ETMA), as well as effects of a harmonic trap within the 
local density approximation (LDA), we showed that, near T c , the local spin susceptibility is 
anomalously suppressed in the trap center due to the formation of preformed singlet Cooper 
pairs. The spatial region where this spin-gap phenomenon occurs becomes wide with de¬ 
creasing the temperature. We also confirmed that the so-called universal thermodynamics 
hold for the spin susceptibility. We pointed out that, using this, we can determine the spin- 
gap temperature T s in a uniform unitary Fermi from the observation of the spatial variation 
of the local spin susceptibility in the trapped case. 

In this paper, we have treated effects of a harmonic trap within LDA, where spatial 
correlations are completely ignored. In addition, the present analyses is restricted to the 
unitarity limit. Improving these issues remains as our future problems. Since a real ultracold 
Fermi gas is always trapped in a harmonic potential, our results would be useful for the 
study of how the spatial inhomogeneity affects thermodynamic properties of this system in 
the BCS-BEC crossover region, as well as how to observe the spin gap temperature T s in a 
unitary Fermi gas. 
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